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Abstract 

The ultra high energy cosmic neutrinos are powerful astrophysical 
probes for both astrophysical mechanisms of particle acceleration and 
fundamental interactions. They open a window into the very distant and 
high-energy Universe that is difficult to access by any human means and 
devices. The possibility of detecting them in large exposure space-based 
apparatus, like JEM-EUSO, is an experimental challenge. In this paper 
we present an estimation of the feasibility of detection of UHE tau neu¬ 
trino by the JEM-EUSO telescope. The interactions of tau-neutrino in 
sea water and Earth’s crust have been investigated. The estimation of 
the propagation length and energy of the outgoing tau-lepton shows that 
if its decay occurs in the atmosphere close enough to the Earth’s surface, 
e.g. below ~ 5km altitude, the cascade is intensive enough and the gener¬ 
ated light can be detected from space. We have evaluated the geometrical 
aperture of the JEM-EUSO detector for the Earth-skimming (horizon¬ 
tal and upward-going) tau-nentrinos by making specific modifications to 
the standard CORSIKA code and developing an interface to the existing 
ESAF (EUSO Simulation and Analysis Framework) software. 


‘Corresponding author: galia@phys.iini-sofia.bg 


1 



1 Introduction 


The investigation of the ultra high energy cosmic rays (UHECRs) and 
neutrinos origin is one of the main scientific objectives of the experiments 
recently searching for high energy cosmic neutrinos such as ANITA [T], 
IceCube [2], the Pierre Auger Observatory [3], RICE [3] and the forthcom¬ 
ing ARA [5], ARIANNA [5], JEM-EUSO [7] and LUNASKA [g]. Since 
the attempts of correlating the arrival directions of UHECRs with known 
active galactic nuclei (AGNs) yield no convincing results [51ll0| other high 
energy sources to which the existence of the UHECRs could be connected 
have been suggested. These include hypernovae (HNe), galactic shocks, 
gamma-ray bursts, etc. A possible clue for the sources of interest would 
be the detection of UHE neutrinos. 

The UHE neutrinos are expected to be born as muon or electron neu¬ 
trinos. Due to vacuum oscillations, however, the neutrino flux at the 
Earth is expected to be almost equally distributed among the three neu¬ 
trino flavours. Neutrinos can interact in a given medium either through 
charged current (CC) or neutral current (NC) channels: viN —>■ X and 
viN —>■ viX, where I = e, fi, t is the outgoing lepton and X is the hadron 
final state. The CC interactions are more important for the space ob¬ 
servations because in the NC interactions the secondary neutrino takes 
away most of the incident energy. Thus, the shower produced by the low 
energy hadronic component resultant from the NC interaction can not 
be observed by a detector configuration such as JEM-EUSO, because the 
shower energy is typically lower than the energy threshold of the detec¬ 
tor except for the very high-energy upcoming neutrinos (of the order of 
lO^^eU). 

Approximately two out of three tau-leptons decay hadronically m, 
Table [T] Erom the leptonic channels only r e u^Ut is relevant for 
tau-lepton detection through initiation of extensive air shower. Having a 
muon instead of electron in the decay would make the tau-lepton escap¬ 
ing detection due to the large decay length of the muon at the energies 
considered thus the shower produced cannot be detected from space. 


Decay channel 

Secondaries 

Branching ratio [%] 

T~ —>■ 7 r~ Vt 

7T~ 

11.8 

T~ TT~ TX^Vr 

n~, 7 r° —>■ 27 

25.8 

T~ —>■ n~ 2n^i'T 

n~ , 2n^ —> 47 

10.79 

T~ — >■ TT~ TX~'k'^Vt 

2n~, 

10.0 

T~ —>■ 7 r~ 'R~ 

2n~, TT^, 7 r° — >■ 27 

5.18 

T~ —>■ n~ R^-K^TV^I't 

27 r“, Stt® — >■ 67 

1.23 

Total hadronic 


64.8 

T~ — >■ e~ Ve.V-r 

e~ 

17.8 

T~ — >■ fl~Uel'T 

h" 

17.4 

Total leptonic 


35.2 


Table 1: Main tau lepton decay channels. 

A cosmic tau-neutrino could interact in the Earth’s atmosphere, sea 
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water or Earth’s crust. It has been shown that the detection of Earth- 
skimming tau neutrinos in the E > lO^^eV range is possible due to the 
fact that the tau-lepton originating in the neutrino CC interaction can 
decay close to the Earth’s surface and the extensive air shower produced 
is intensive enough to be observed by both fluorescence and Cherenkov 
detectors misi m US]. In the atmosphere, the neutrino interaction 
probability increases with the zenith angle. The Earth’s curvature and 
the density of the different atmospheric layers are the key parameters 
that influence the neutrino induced extensive air shower development. A 
very high energy tau-neutrino could interact deeply in the atmosphere 
and, instead of one-peak signature, to initiate the remarkable double¬ 
bang event that could be observed by space based telescopes with the 
JEM-EUSO design parameters. 

In this work we study, by simulations, another possibility, first pointed 
out in na, namely Interactions of quasi-horizontal 85'^ < 6 < 90'^ and 
upward going 90° < 9 < 95° tau-neutrinos in higher (than air) den sity 
media such as Earth’s crust and see water. In the simulations we track out 
the decays of tau-leptons emerging from ocean or Earth’s surface in order 
to estimate the feasibility of detection of such events by the JEM-EUSO 
telescope. 

1.1 JEM-EUSO telescope 

JEM-EUSO, Extreme Universe Space Observatory (EUSO) attached to 
the Japanese Experiment Module (JEM) on board the International Space 
Station (ISS), is a future mission devoted to investigation of extreme en¬ 
ergetic cosmic rays. The JEM-EUSO apparatus is a remote-sensing space 
instrument |16| . Fig[T]that will orbit the Earth every « 90 minutes on the 
ISS at an altitude Hiss ~ 400 km. It consists of photon collecting optics 
|17| and focal surface (FS) detector |18| together with its electronics m- 
The telescope optics is composed of three double-sided curved circular 
Fresnel lenses with 2.65 m maximal diameter. When high energy cosmic 
rays interact with the atomic nuclei of air molecules they initiate exten¬ 
sive air showers. JEM-EUSO will detect the light from isotropic nitrogen 
fluorescence excited by the extensive air showers and Cherenkov radiation 
reflected from the earth surface or dense clouds. When the ultraviolet 
light hits the optics it is focused onto the focal surface which consists of 
137 Photo-Detector Modules (PDMs). Each PDM consists of a set of el¬ 
ementary cells containing arrays of Multi-Anode Photo Multiplier Tubes 
with 64 pixels with a spatial resolution of 0.074° each. Data ftom each 
pixel, is acquired within 2.5 fis time window when the trigger logic (de¬ 
fined according to the hardware requirements) is activated. The intensity 
of the observed light depends on the transmittance of the atmosphere, 
the cloud coverage and the height of the clouds. For precise measure¬ 
ments JEM-EUSO is equipped with an Atmospheric Monitoring system 
m that will characterize the atmospheric conditions and thus determine 
the effective observation aperture with high accuracy. The atmospheric 
monitoring system consists of an infrared camera and a light detection and 
ranging system. The observation area of the earth surface is determined 
by the projection of the field of view (FoV) of the optics and the area of 
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the focal surface. The dimensions of the FoV are 64° and ~ 45° on the 
major and minor axes, respectively. For these axes, the projected lengths 
on earth surface are 500 km and ~ 330 km, respectively for altitude of 
the ISS orbit Hiss = 400 km. The effective solid angle QfoV is ~ 0.85 
sterad. 

For the baseline layout of 137 PDMs, the observation area of the ap¬ 
paratus, for nadir mode, is a function of Hiss expressed by: 




Figure 1: Observational principle of the detector. H is the orbital height, 7 m 
the FoV aperture semi-angle and atm is the tilt angle between the optical axis 
and the nadir. 



( 1 ) 


For tilt angles a < 40°, Aobs{a) = 71^^“'^*’^^(cos(q))“*’, where, for the 
altitude of interest b is in the range 3.2 < b < 3.4. For a ~ 40° — 50°, the 
detector is sensitive to the area over the local horizon, and Aoba saturates 
above a ~ 60°. 

2 Simulation of neutrino interactions in 
water and earth crust 

To study the capability of the JEM-EUSO detector to detect UHE neu¬ 
trino events we have adapted the version 7.401 of CORSIKA code EH 
with new options for neutrino interactions in water and earth media. For 
this we have added one layer below the lowest atmospheric one. It cor¬ 
responds to the last or external layer of the Earth (Earth’s crust) with 
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density profile chosen according to the Preliminary Earth Model |22| or 
to a layer with sea water density p-u, = 1.025(//cm^. The user is free to 
choose the thickness of the layer and therefore to define the interaction 
volume. The observation of tau-neutrinos depends on the incident fluxes 
(for which a large number of models exists), the cross section and on the 
detector capability. The major uncertainty on the cross section at the 
energies we consider comes from the unknown behaviour of the parton 
distribution functions (PDFs) at very small values of the parton momen¬ 
tum fraction x. For the simulation of the charged and neutral current 
interactions we have used the deep inelastic neutrino-nucleon cross sec¬ 
tions |23| . see Fig. [5] and extrapolated to very small x parton distribution 
functions CTIO |23i 1241125| . The high energy cross-section is evaluated 
specifically for experiments devoted to the UHE neutrino measurements. 
The input PDFs considered treat heavy quarks by using general-mass- 
variable-flavour number schemes. 

For primary neutrinos CORSIKA calls the HERWIG Monte Carlo gen¬ 
erator |26| to produce the first interaction and then tracks the products 
of this interaction. The call to HERWIG has been adapted to the cases of 
water and earth taking into account the dependence of the interactions on 
the inelasticity parameter y — 1 — In absence of reliable predictions 
for the shape of the UHE neutrino spectrum, we consider power law fluxes 
with 7 = 2 to study how the propagation of tau is influenced 
by its energy losses. 






Figure 2: Neutrino DIS cross-section for CC and NC scattering as predicted by 
the HERAPDF1.5. and total relative PDF uncertainties including experimental 
errors. [23] 
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Ntau/Nnutau 


2.1 Tau lepton propagation and energy losses 

The average lepton energy loss per distance X [g/cm^] traveled is a sum 
of nearly constant ionization energy loss and weakly energy dependent 
radiative energy loss through bremsstrahlung, pair production and photo- 
nuclear scattering [27|[^|^|^[3l]. For the case of earth medium and 
energy region E > 10^^ eV the dominating photo-nuclear process is re¬ 
sponsible for the most uncertain contribution. In our code the energy 
losses of the charged leptons in water are computed by a modified stop¬ 
ping power formula to account for the density effect. The cross sections 
for tau-lepton interactions with mater we use are taken from |2711291130| . 
In these cross sections the nuclear structure of the target (elastic and in¬ 
elastic form factors), screening (shadowing) effects are taken into account. 
Following the respective CORSIKA procedures we treat the tau propaga¬ 
tion stochastically, with a series of electromagnetic interactions occurring 
according to the cross section and decay probability. We have simulated 
10® tau-neutrino events in both media with energies in the range 
10^® — 10^^ eF. Our simulations show that the tau-lepton, produced by 
CC Vt interactions under the surface of the Earth or the ocean emerges 
in the atmosphere with energies shown in Fig. [S] 




Figure 3: Energy distribution of tau-leptons emerging from earth (left) and 
water (right) surfaces normalized to the number of simulated events. 


The distance R (propagation length) that tau travels after its pro¬ 
duction in earth or water medium as a function of the incident neutrino 
energy is shown in Fig. |4] For energies above the JEM-EUSO threshold 
of 3 X lO^^eV its propagation length ranges from ~10 to ~15 km. Our 
results are in a good agreement with the results of [M131 [33] despite the 
difference in methods for calculation of R that have been used in those 
works. 

The number of showers produced, for tau-neutrino energies above E = 
3 X lO^^eF, is shown in Table [2| As one can see the number of showers 
decreases with increasing of the primary neutrino energy. For a shower 
that provides a detectable signal on the focal surface of the detector the 
energy of its charged component has to be larger than Echargsd ~ 2.5 x 
10'^®eF. To test the JEM-EUSO sensitivity to the skimming events we 
have analyzed all the showers with energy of the charged component larger 
than 2 x 10^^eV. The results indicate that if the decay of the emerging tau- 
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Figure 4: Distance R (propagation length) that the tau-lepton travels in earth 
(blue line) and water (red line) medium after its production as function of 
incident neutrino energy. 


xlO""' [eV] 

3-5 

5-7 

7-9 

9-11 

11 - 13 

interaction in earth 

0z = [85^ - 90“] 
e, = [90° - 95°] 

2230 [895] 
3864 [741] 

1012 [689] 
1484 [663] 

629 [99] 
856 [102] 

493 [49] 
615 [52] 

415 [0] 
507 [1] 

interaction in water 

e, = [85" - 90"] 
e, = [90° - 95°] 

2291 [1017] 
3786 [984] 

1315 [689] 
1341 [663] 

714 [102] 
805 [106] 

513 [54] 
621 [58] 

442 [1] 
512 [2] 


Table 2: Number of showers produced in the atmosphere for tau-neutrino ener¬ 
gies above E = 3 x lO^^eV. The number of the initial interactions is 10® and of 
the emerged tau leptons is ^ 25 x 10®. In brackets: A selection of the showers 
generated from tau decays occurred below H = 6 km and Echarged > 2 x 10®® 
has been applied. 


lepton occurs in the atmosphere close enough to the Earth’s surface, e.g. 
below ~ 5km altitude, the cascade is intensive enough and the generated 
light can be detected from space. The number of showers passing the 
selection criterion is presented in square brackets. Table [2l 

3 Estimation of the geometrical aperture 

The estimation of the aperture has been done by means of ESAF |34| 
software and a special interface of this code to the CORSIKA package. 
Both software parts are written in CH—h and the output is visualized by the 
ROOT package |35| . In our analysis we assume clear sky condition with 
average background. By rescaling measurements of the Tatiana satellite 
|36| to the ISS orbit |37l I38| the average background flux is estimated to 
be ~ 500 photons.m~^.steradr^.ns~^. 

The acceptance of the detector is: 
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Acc= I I dSdQe{l^, 6, E, k)cos9, 

Js Jn 

where 9, E, k) is the detection efficiency depending on the EAS posi¬ 
tion 1^ , zenith angle 9, energy E and all possible causes of a lack of detec¬ 
tion that we denote by k. We determine the geometrical aperture through 
the probability of satisfying the second level trigger (Linear Track Trigger) 
m as obtained from our simulation. The trigger follows the movement 
of the shower spot in the focal surface over a specific time window to 
distinguish the event pattern from the background. To take into account 
also the showers crossing the field of view with a core location out of the 
observation area, we have injected a number of showers, Ninj, initiated 
by tau-leptons with energies shown in Fig. [3] and emerging angles close 
to 90'^ in an effective area Sinj bigger than the observation area of the 
telescope. 

For the present JEM-EUSO baseline configuration the geometrical 
aperture of JEM-EUSO as a function ot the primary energy is: 

(3) 

where the ratio gives the probability of a trigger signal to be 

created when a shower location V' is inside the detector field-of-view. Uq 
is the effective solid angle of the detector for the zenith angle range of in¬ 
terest. Ninj is the injected number of showers initiated in the atmosphere 
by the tau-leptons emerging from the water or earth medium and Ntrigg 
is the number of the events selected by the trigger logic mentioned above. 

Fig. [5] shows the geometrical aperture as a function of the primary 
energy for the case of horizontal and upward going neutrinos. 
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Figure 5: Geometrical aperture of the JEM-EUSO for horizontal and upward 
going neutrinos interacted in earth (left) and water (right) medium. 

For both media the geometrical aperture increases with the neutrino 
energy. In case of interactions in the air the dependence of the geometrical 
aperture on the neutrino energy is approximately the same but it follows 
the trend until higher energies. It is due to the smaller energy losses of 
tau and the Earth’s curvature; both influence the propagation length of 






the lepton and the longitudinal particle and energy distributions of the 
showers produced. Estimation of the detector acceptance and expected 
event rate for tau neutrino interactions in the three media is subject of a 
forthcoming paper. 


4 Conclusion 

A tool for simulation of UHE cosmic tau-neutrino interactions in water 
and Earth’s crust has been developed. The propagation length of the pro¬ 
duced tau-leptons in both media as well as their energy distributions are 
presented. The simulation of the detector response to the fluorescent and 
Cherenkov components of the extensive air showers has been performed 
by the ESAF package with the interface to the CORSIKA software that 
we have developed. The results shown in Fig. 3 indicate that the trigger 
efficiency of the JEM-EUSO telescope increases up to ~ 100% at ener¬ 
gies > 8 X lO^^eE. The geometrical aperture of JEM-EUSO detector 
for Earth-skimming tau-neutrino events has been estimated for clear sky 
condition, nadir mode and average background. The detected events at 
energies close to the detector threshold will make possible a cross-check 
based on the obtained energy spectrum with the ground-based experi¬ 
ments. The detailed simulation analysis of the event rate per year and 
full operational time for tau-neutrino interactions in atmosphere, sea wa¬ 
ter and Earth’s crust is in progress. 
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